The discovery of an in situ Neanderthal remain in the Bawa Yawan Rockshelter, West-Central Zagros Mountains, Kermanshah

Neanderthal extinction has been a matter of debate for many years. New discoveries, better chronologies and genomic evidence have done much to clarify some of the issues. This evidence suggests that Neanderthals became extinct around 40,000–37,000 years before present (BP), after a period of coexistence with Homo sapiens of several millennia, involving biological and cultural interactions between the two groups. However, the bulk of this evidence relates to Western Eurasia, and recent work in Central Asia and Siberia has shown that there is considerable local variation. Southwestern Asia, despite having a number of significant Neanderthal remains, has not played a major part in the debate over extinction. Here we report a Neanderthal deciduous canine from the site of Bawa Yawan in the West-Central Zagros Mountains of Iran. The tooth is associated with Zagros Mousterian lithics, and its context is preliminary dated to between ~43,600 and ~41,500 years ago.


Introduction
Neanderthals were a very successful hominin lineage that existed for several hundred thousand years, and their extinction remains one of the most persistent questions in palaeoanthropology. With such a vast geographical range, which at times extended from westernmost Europe and across a very large area of Asia [1,2] identifying the last surviving Neanderthal populations is critical for interpreting the mechanisms behind their demographic decline, the conditions that enhanced or reduced their resilience, and the geography and timing of their last interactions with modern humans. Unsurprisingly, one of the most controversial aspects of a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 this debate is the chronology of their disappearance. Aspects of that chronology are relatively well established, in Europe, where Neanderthals appear to have a continuous occupation since the Middle Pleistocene, their last occurrence has been timed to~40,000-37,000 cal. BP [1][2][3][4]. A number of sites testify to the presence of Homo sapiens in Europe from at least 46,000 BP-Bacho Kiro [1] and Grotta del Cavallo [5,6]. Therefore, the chronology of Neanderthal extinction across most, if not all of Europe, implies a minimum of 4,000 years of contemporaneity with Homo sapiens groups [1,2], during which both biological [7] and cultural [8] interactions took place. At the other extreme of their geographical range, the Neanderthal occupation of parts of Siberia reveals a very different pattern. Palaeoanthropological and ancient genomic evidence from Denisova, Chagyrskaya and Okladnikov caves suggests repeated west-to-east dispersals between 150,000 and 120,000 years [9][10][11] which gave rise to at least two temporary Neanderthal populations that became demographically vulnerable, characterised by episodes of interbreeding with local Denisovans and small population sizes [11][12][13].
However, there is one region that presents a particular challenge because of the sheer complexity of both its prehistoric record and its geographic setting-southwest Asia. Bordered by the Mediterranean Sea in the west and mountain ranges in the north and east, the core of the region consists of deserts and semi-deserts that, even during humid climate intervals, would have constrained resources along waterways [14]. Thus, the prehistoric record of southwest Asia is concentrated along an arc extending from the Levant in the west to the Caucasus in the north and the Zagros Mountains to the east. These areas-Levant, Caucasus and Zagros-in turn, offered hominins different challenges and opportunities, and this is reflected in their different records.
Biogeographically, the Levant oscillated between African and Eurasian biomes, acting as a temporary corridor for dispersing African populations [15]. It is now known that this corridor was used by Homo sapiens at least twice-between 177,000-194,000 years, as evidenced at the site of Misliya [16] and between~120,000 and 90,000 years, as shown at the sites of Skhul and Qafzeh [17], before the area was permanently occupied by Homo sapiens around 55,000 years [18]. These intermittent early modern human occupations of the Levant were interspersed by the use of the area by Neanderthals [19], whose local extinction is represented by the last Middle Palaeolithic (MP) levels at Kebara at 48,000-49,000 BP [20]. However, the presence of Homo sapiens at Manot Cave~55,000 years [18] suggests a relatively long period of overlap between the two groups in the area.
In contrast, the Caucasus, a 1200 km mountain range between the Black and Caspian Seas, was an important barrier to animal and hominin dispersals north, as shown by the genetic differentiation of several taxa into southern and northern Caucasus populations, and also reflected in differences in lithic technology [21,22]. Southern Caucasus MP assemblages appear to have been part of a Levantine network [21] that extended to the Zagros during the terminal MP phases [23]. However, in contrast to the Levant, Neanderthals occupied the Caucasus for nearly 10,000 years longer, with an age for the latest Neanderthal and MP industries similar to that of European Neanderthals, at~40,000 cal. BP [22]. Yet, differently from both the Levantine and European records, Homo sapiens occupation of the Caucasus does not appear to occur prior to the local extinction of Neanderthals, thus suggesting no local contemporaneity of the two groups [22]. Ultimately, another picture of Neanderthal-Homo sapiens landscape interactions emerged from the Zagros Mountains which stretch over 1600 km in a northwestsoutheast direction and reaching 4400 m above sea level (asl) in Iran. Moreover, due to the steep and rugged topographic conditions, the Zagros embrace diverse ecological ecotones. They represent a formidable geographical barrier, but also a diverse set of habitats and ecozones [24] that were exploited by hominins since the Lower Palaeolithic [25]. MP sites are relatively numerous, although only three have yielded Neanderthal fossils. The best-known of these is Shanidar, where the remains of ten Neanderthals were discovered [26][27][28]. The Shanidar Neanderthals have played a major role in the discussions about the complexity of Neanderthal behaviour, both for the presence of older-aged individuals in the fossil assemblage and the controversial claims of a 'flower burial'. Yet, despite their importance, their timing and contemporaneity have not been finely resolved; a broad age between 70,000 and 45,000 years has usually been reported, although recent excavations have suggested an age range of 55,000-45,000 years [28]. Approximately 350 km to the southeast, among a number of MP and UP sites, Wezmeh and Bisetun Caves in the Kermanshah region (Fig 1) have also produced Neanderthal remains. Wezmeh Cave has no artefactual evidence for human occupation in the Late Pleistocene, and the presence of a human premolar tooth (Wezmeh 1), recently shown to have Neanderthal affinities [29], was interpreted as the result of carnivore hunting/scavenging activities. The age estimation for this premolar is very poor, bracketed between 70,000 and 11,000 years based on the uranium-series analyses of the fauna by alpha spectrometry [29].
Nevertheless, the transitional period from Yafteh, Ghār-e Boof and Shanidar cave is bracketed between 45,100 and 40,350 cal. BP, with a 68.2% probability [30,31]. Together with the approximate youngest age of the Neanderthal remains from Shanidar at 45,000 BP, this would suggest a pattern similar to the southern Caucasus, where modern humans succeed, rather than overlap in time with the local Neanderthals. Moreover, a probable, but undated, Neanderthal right proximal radius diaphysis was found in Bisetun Cave in association with typical Zagros Mousterian lithic artifacts [31,32].
Despite these few Neanderthal fossil hominin remains, there is a widely held assumption that the MP assemblages in the Zagros Mountains were produced by Neanderthals, and that the local Middle to Upper Palaeolithic transition represents the arrival in the area of Homo sapiens [33]. Here we report the discovery of a new Neanderthal remain at Bawa Yawan rockshelter (BY1) in the Central Zagros Mountains confirming that Neanderthals manufactured Middle Palaeolithic artifacts. The new chronology associate with the site in question is showing a surprising continuity of Neanderthals presence in the area long after the first Upper Palaeolithic Homo sapiens, however, more radiocarbon dates are on the way to strengthen our conclusions.

Materials and methods
All necessary permits for this research, including excavations and materials study to Saman Heydari-Guran, were issued by the Iranian Center for Archaeological Research (Permit number 963141/34/3952). The archaeological materials presented in this research are stored and available for the study with authorization from the Iranian Cultural Heritage and Tourism Organization (ICHTO) in Kermanshah Province.
The results for this study are based on the three excavation seasons between 2016-2018 in the archaeological site of Bawa Yawan (34˚38' 23.70"N, 46˚55' 48.36"E, 1300 m asl), locates in Kermanshah Region, Iran.

μCT images of BY1
"High-resolution μCT images of the lower deciduous canine BY1 were obtained with the following scan parameters: 125 kilovoltage (kVp), 120 μA, and 0.015 mm voxel size. The μCT volume was segmented using Avizo 9.2 software (Thermo Fisher Scientific, Waltham, Massachusetts, US), and the 3D models of the dental tissues (i.e., enamel, dentine and the preserved portion of the pulp chamber) were refined in Geomagic Design X (3D Systems Software, Rock Hill, South Carolina, US) to optimize the triangles and create fully closed surfaces.
The morphological description of BY1 was performed according to standards outlined by the (ASUDAS [34]) Although ASUDAS has been devised for modern human permanent dentition, we applied protocol to fossil deciduous tooth as already done in other publications [e.g., [35][36][37] because it permits a more precise and accurate comparison at each degree of

PLOS ONE
Discovery of an in situ Neanderthal remain in Bawa Yawan development. The wear stage of the occlusal surface was assessed based on [38], while the age of root resorption was estimated according to Moorrees et al [39,40]. The crown diameters of BY1 were virtually measured in Geomagic Design X, orienting the tooth with the best-fit plane computed at the cervical line parallel to the xy-plane of the Cartesian coordinate system, and with the lingual aspect parallel to the x-axis. The size of the bounding box enclosing the crown was used to collect mesiodistal (MD) and buccolingual (BL) diameters.
Comparative data for the BL diameter was collected from the literature [40][41][42]. The computation of the 3D average enamel thickness (AET) and relative enamel thickness (RET) index follow [41]. A set of Neanderthals (NEA; n = 6), early H. sapiens (EHS; n = 2), and recent H. sapiens (RHS; n = 3) lower deciduous canines at different wear stages was acquired using a Skyscan scanner with voxel size ranging from 13 μm to 74 μm) at the Department of Human Evolution of the Max Planck Institute for Evolutionary Anthropology. As for BY1, dental tissues were segmented in Avizo 9.2 software, and the digital models were refined in Geomagic Design X. The 3D AET index (in mm) is the enamel volume divided by the underlying EDJ surface, while the 3D RET index (scale-free) is the 3D AET divided by the cube root of the crown dentine and pulp volumes.

Radiocarbon dating
Thirteen charcoal samples collected from GH 5, 4, 3, and 2 were submitted for 14 C dating to the Klaus-Tsichira-AMS facility of the Curt-Engelhorn Centre in Mannheim, Germany. Prior to dating, the samples were pretreated with the ABOX method and dated [43] in the radiocarbon lab. For samples that yielded finite conventional ages, calendar-year ages (and the corresponding 68.2% and 95.4% confidence intervals) were estimated using the IntCal20 calibration dataset [44].

Nawdarwan Valley, Bawa Yawan rockshelter and stratigraphy
Located in the Nawdarwan Valley (34˚38' 23.70"N, 46˚55' 48.36"E, 1300 m asl), the Bawa Yawan (BY) rockshelter was discovered in 2009-10, during surveys for Palaeolithic sites in the Kermanshah region of the Central Zagros [33]. The Nawdarwan Valley forms a natural 17 km long corridor around the Razawar River that connects the two fertile plains of Kermanshah and Kamyaran in a N-S direction (Fig 1). The valley is rich in MP and UP rockshelters, with more than 50 sites identified during the surveys [33]. One of these, the Bawa Yawan rockshelter, is set in 50 m high limestone cliffs, visible from a far distance, that stand at the edge of a flat plain to the southwest (Fig 2A and 2B). A karstic spring, forming a natural pond, is located 50 m from the site. The archaeological deposits are accumulated on the southern slope of the cliff, occurring on a flat occupational area of ca. 300 m 2 , which is 10 m higher than the plain surface today. Several large boulders lie at the eastern corner of the rockshelter and on the talus slope, protecting the archaeological deposits from erosion ( Fig 3A). Excavations of the Bawa Yawan rockshelter took place during three field seasons between 2016-2018. These focused on two areas, 10 m apart, both located close to the rockshelter wall, named the Western (14 m 2 ) and Eastern (6 m 2 ) trenches ( Fig 3A).
The Bawa Yawan deposits are divided into 5 geological horizons (GH 1-5) based on colour and sediment texture ( Fig 3C). These have a gradient of about 10˚from the back wall towards the outside, but are relatively level in east and west directions (Fig 3A and 3B). The deposits are mostly composed of clastic materials, and fine sediments washed into the site through joints and bedding planes from the back-wall cliff. In some squares (H31, I31) the excavation reached a maximum depth of 4.5 m from the surface without reaching bedrock (Fig 3B and 3C). Some horizons, particularly GH2 and GH5, are characterised by a significant number of irregular stones collapsed from the cliff, varying in size from large boulders to small pebbles. GH3, with an average thickness of 80 cm, is different from the other geological horizons in being harder and having a higher carbonate content, protecting the underlying deposits from erosion.
• GH5: a dark reddish-brown (dry 5YR3/4 and moist 5YR2/4 very dark reddish-brown) and only exposed in five squares at Western trench and with 150 cm includes the lowermost geological layer and the deepest part of the excavation (Fig 3B and 3C). The sediments' colour becomes darker containing a lot of iron oxide particles, specifically in the lower part of the section. The deposits range from sandy/silty granules up to angular limestone of different sizes from small pebbles to boulders coming from the back wall of the shelter. The contact between GH4 and GH5 is gradual. The sediments become gradually reddish-dark brown downwards, specifically in the lowermost part of the section. It is the richest layer throughout the stratigraphy and is associated with numerous bone fragments which are partly burnt, very small charcoal fragments, and MP lithic artifacts of the Zagros Mousterian tradition. This GH yielded a good amount of fauna including Bovidae, Bovidae antilopinae and Equidae. So far, bedrock has not been reached during the three excavation seasons and at a depth of 4.70 m (Fig 3).
• GH4: This GH with about 65 cm thickness and a light reddish-brown colour (dry 5YR4/6 reddish-brown and moist is 5YR7/3 dull orange) stands in sharp contrast with the underlying horizon (GH5). The amount of angular cobbles and boulders is more than the upper GHs. This GH faces a dramatic decrease in cultural materials, indicating a sudden abandonment of the site. The few lithic artifacts recovered in this layer belong to MP period (Fig 3).
• GH3: Around 95 cm thick, the GH3, changes relatively gradually from GH4 (dry 5YR7/3 dull reddish-brown and moist 5YR ¾ dark reddish-brown). It contains light colouartr compacted highly calcareous deposits with concentration of eboulis and angular pebbles. The number of lithic artifacts drops dramatically in this GH. Most of them are typical MP artifacts, some of which indicate Levallois technology. Very few UP artifacts have been also observed but only at the uppermost of GH3. The faunal remains are scattered and very fragmented as well (Fig 3B and 3C).
• GH2: This layer was observed in both Western and Eastern trenches. Up to 80-95 cm thick, GH2 is in general composed of reddish-brown deposits (dry 2.5YR5/4, moist 2.5YR2.5/ 3-dark reddish-brown) mixture of sand-silt-clay grain sediments associated with cobbles and medium-sized pebbles. This layer yielded archaeological materials from three different periods of MP, UP and Epipaleolithic (Fig 3B and 3C). GH2 yielded relatively abundant fauna, but highly fragmented including Bovidae, caprinae, and small birds.
• GH1: The uppermost geological layer is dull reddish-brown (dry 5Y6/4 dull orange, moist 2.5YR4/4, dullreddish-brown) and is around 30 cm thick, containing the recent soil associated with a few mixed modern objects and it is separated by an erosional unconformity to the underlying layer (GH2).

BY1 Neanderthal tooth.
At a depth of around 2.50 m from the surface in the upper part of GH5, where the density of MP cultural materials gradually reduced, an in situ hominin tooth (hereafter called BY1) was recovered in association with fauna and Zagros Mousterian lithic artifacts (Fig 3B and 3C). The tooth is an exfoliated lower left deciduous canine (Ldc1) consisting of a relatively well-preserved crown and about one-fourth of the root (Fig 4). Neither caries nor enamel hypoplasia is visible. The enamel shows several longitudinal fractures from the cervix to the incisal surface that affect the underlying dentine. The incisal surface is worn obliquely, mesially to distally, exposing a large area of dentine up to wear stage 4 [38]. The buccal surface exhibits a mesiodistal convexity with its maximum at the mesial aspect (Arizona State University Dental Anthropology System (ASUDAS) grade 4). Accordingly, in the incisal view, the crown appears asymmetrical, which is further emphasized by the distal projection of a moderate lingual cervical eminence. The lingual surface is concave, bordered by moderately expressed distal and mesial marginal ridges (as also clearly shown in the EDJ, Fig 4). These ridges merge at the cervical eminence giving a semi-shovel shaped aspect to the crown (ASUDAS grade 4).
The moderate buccal bulging of the crown and the slight flaring of the mesial and distal sides from the cervix are more consistent with a lower than an upper deciduous canine. Interproximal wear facets are visible on both mesial and distal sides, the distal side being smaller (length: 1.20 mm; height: 1.73 mm) than the mesial one (length: 1.49 mm; height: 2.59 mm). The preserved root, slightly more elongated labially (mid-labial height: 2.95 mm) than lingually, (mid-lingual height: 2.31 mm) is resorbed (stage Res3/4 [45]), suggesting an age at exfoliation of approximately six years on the basis of recent human standards [38]. The crown has a mesiodistal (MD) diameter of 6.87 mm (minimum estimation due to wear) and a buccolingual (BL) diameter of 6.50 mm, which is close to UP Homo sapiens and, particularly, Neanderthal mean values (Table 1). At the cervix, the MD and BL diameters are, respectively, 5.32 mm and 5.88 mm. Despite the reduced comparative sample and the intrinsic limits due to the wear stage of the teeth, unworn Neanderthal lower deciduous canines show the average and relative enamel thickness indices (AET and RET respectively) lower than the values observed for worn early and recent Homo sapiens specimens ( Table 2). The Z-scores computed for the AET and RET values of the BY1 specimen fall very close to the mean values of Neanderthals at similar wear stage (i.e., Krapina 51, Tagliente 4), further supporting its Neanderthal affinity.
Cultural materials and artefact technology and typology. The technological and typological characteristics of the Bawa Yawan artifacts are based on detailed studies of the lithic assemblages from GHs 5 to 2. Three techno-complexes of MP Mousterian, UP "LaK" [46] and the Epipalaeolithic from Bawa Yawan have been compared with their counterparts throughout the Zagros Mountains (Figs 5-7). We use attribute analysis and chaînes opératoire approaches to compare specifically the assemblages from the West-Central Zagros sites including Warwasi and Shanidar. Key attributes include core preparation and reduction strategies, degree of core reduction, modes of flaking, platform preparations, and the retouching and re-sharpening of the tools. We recorded the technological attributes of each artefact individually and in relation to the reduction sequence in order to reconstruct the chaînes opératoires in each archaeological horizon.  The excavations revealed three archaeological horizons throughout the stratigraphy, resulting in a total of 10,411 lithic artifacts (Tables 3 and 4); those artifacts that were in primary depositional context were concentrated in two dense bands in GH5 (MP) and GH2 (MP, UP and Epipalaeolithic) (Fig 3B). MP materials were recovered from GH3 and GH4, although at lower concentration levels. At the top of the sequence in the Eastern trench, close to the rockshelter wall, there is a rich Epipalaeolithic industry, similar to that reported at the Warwasi rockshelter [47], characterised mainly by geometric microliths and microburins. Following the Epipalaeolithic, the first half of GH2 in both the Eastern and Western trenches preserves UP lithics typical of the West-Central Zagros, especially in Warwasi [46]. The lithics are characterised by laminar technology. The tools mainly dominated by different kinds of twisted and retouched bladelets and end scrapers. At the bottom of GH2, and throughout GH3 to GH5, in both the Eastern and Western trenches, lithics typical of the MP Zagros techno-complex, namely the Zagros Mousterian [48], were uncovered (Tables 3 and 4

, Figs 5 and 6).
Chronology. Thirteen charcoal samples were collected from different parts of the Bawa Yawan stratigraphy: (a) six charcoal samples were obtained from the below and above the hominin tooth in GH4 and GH5; (b) six samples from squares I29 and G27 from the beginning and middle of the UP occupation in GH2 and GH3; (Fig 8) (c) and one sample from the Epipalaeolithic layer at the upper part of GH2 in the Eastern Trench.
Eight of the twelve samples were successfully dated, while one returned a minimum date (Tables 5 and 6). All the 14 C results are calibrated using the new IntCal20 [43] within the OxCal 4.4 program [43,49] and are discussed here at 68.2% probability. The youngest occupation of the site, corresponding to the Epipalaeolithic level, ranging between 13,400 and 13,300 cal. BP (Fig 9) Only two of the six charcoal samples from the UP GH2 were dated successfully (both from square G27: MAMS-41723 and MAMS-41725). The two ages obtained are consistent with each other, yielding calibrated age ranges of 34,700-34,400 and 34,600-34,400 cal. BP, respectively ( Fig 3B and Table 5). Five of the six charcoal samples collected from MP occupation layers (GH4 and GH5) produced reliable radiocarbon ages; a fifth sample produced a minimum age. (Table 7 and Fig 8). The ages obtained are consistent with their respective depths within the stratigraphy. The sample from the middle part of GH4 (MAMS-44523), 30 cm above the level of the BY1 in the same square (I31), provided an age of 40,300-39,400 cal. BP. Another sample from square I31 (MAMS-44525), at 9 cm above BY1, provided an age between 42,100 and 41,800 cal. BP, and a charcoal sample at a level just 2 cm above BY1 in square H31 (Beta-465001) provides the closest age estimation of the specimen, ranging between 43,000 and Table 2 (Table 7).

. Values of the components of enamel thickness computed for Neanderthal (NEA), early Homo sapiens (EHS) and recent Homo sapiens (RHS) lower deciduous canines at different wear stages (ws).
We performed a Bayesian model using the new IntCal20 in OxCal 4.4 and the t-type outlier analysis (prior probability set at 5%), to detect if there is any problem between the 14 C ages and their relative stratigraphic position at the site (Fig 9 and Table 6). The model, which was run

PLOS ONE
Discovery of an in situ Neanderthal remain in Bawa Yawan using only the dates from GH5 until the GH2, confirm the integrity of the chronology with an Agreement index of 107%, well above 60% (Table 7 and Fig 8). The layer GH5 start at 44,600-43,100 ending in a transitional boundary with the GH4 at 42,000-40,500 cal. BP. In reality, one 14 C age at the bottom of the sequence, result in an age older than 44,200 14 C BP, hence it could be older than the ages provided by the Bayesian model. Using the 'date' command in OxCal we estimated the duration of each phase at the site, and knowing the exact stratigraphic depth of the Neanderthal tooth in relation with all the others samples, an indirect date of BY1 could have been determined (see Method section). The total duration of GH5 ranges between 43,600 and 41,500 cal. BP, the duration of GH4 is between 41,100 and 39,000 cal. BP, with the upper part (GH2) ranging between 34,900 and 34,200 cal. BP (Table 7).
These dates indicate that the Bawa Yawan rockshelter contains a sequence from the MP to the UP in both the Western and Eastern trenches, as well as an Epipalaeolithic occupation level above those in the Eastern trench, which was not modelled. According to the Bayesian age modelling, GH4 and GH3 are broadly simultaneous with a substantially cold period in the Northern Hemisphere-the so-called Heinrich Stadial 4 (Fig 8), which might be responsible for the reduced cultural materials in this GH.
Furthermore, at this stage of research, these ranges show two critical ages for the Neanderthal occupation of the Zagros Mountains-first, the age of the Neanderthal BY1 tooth can be bracketed between 43,600 and 41,500 cal. BP, and second, the possible age of a young Zagros Mousterian in the area around 39,000 cal. BP. It is important to stress that more radiocarbon dates are on the way to help in this interpretation since the young age in GH2 is based only on a single sample.

Discussion
On current evidence, we can identify multiple events of Neanderthal local extinctions-first, a Neanderthal group, unrelated to later local populations, becomes extinct in the Altai during Marine Isotope Stage (MIS) 5 (130,000-74,000 BP) [11,12], second, the extinction of a The lithic artifacts are divided into technological categories based on artificial 10 cm spit. AD: angular debris, SD: small debitage (> 10mm), MD: micro debitage (> 5mm).
https://doi.org/10.1371/journal.pone.0253708.t003   [2,4,22]. Contrasting with the synchronicity of their final disappearance, these late Neanderthal groups appear to be spatially structured and to have different demographic trajectories. The discoveries from the Bawa Yawan rockshelter bring clarity to the position of the Zagros Mountains in this emerging picture. The presence of Neanderthals in the Zagros has been known for decades [26,27,29,32]. However, few sites (Shanidar, Warwasi, Ghar-e-Khar, Kaldar, Ghamari, Gilvaran) preserve both MP and UP industries [25-27, 51, 52], while chronometric ages derive only from the old and new excavations at Shanidar, which point to the use of the cave by Neanderthals between 70,000 and 45,000 BP [26][27][28]. The age between 43,000 and 41,000 cal. BP for BY1 tooth associated with the Zagros Mousterain artifcts in the West-   Yawan together with the information from southern Caucasus [21,23] suggest that caves and rockshelters throughout the entire mountains arch that extends from Georgia to Iran were used by relict population of Neanderthals could survived a short time after they had become extinct across most of Europe. Furthermore, similar to the records of the Levant and Europe, the Zagros early UP assemblages are techno-typologically diverse, showing variable patterns of similarities and differences to early UP industries elsewhere, and suggesting a diversity of modern human groups with varying foraging strategies and ranges along the vast mountain chain [53]. The Bawa Yawan discoveries hypothesizes that the final disappearance of Neanderthals took place in the context of increasing competition with Homo sapiens, the multiple events of localised Neanderthal extinction, together with increasing genomic evidence of comparatively smaller population sizes in marginal areas [11]. The combination of these differences in demographic resilience between the two groups during periods of repeated climatic insults [54,55], coupled with the more sophisticated technology of expanding modern humans (e.g., [56,57] may have been the tipping point that created what appears as a vast spatial extinction event across thousands of kilometres, but was the sum of numerous local independent population processes, among which the disappearance of Neanderthals in the Zagros Mountains was one of the last. Finally, it will always be the case that artefactual data will be far more abundant than fossil evidence. Yet, as recent work in Central Asia and Siberia has shown [12], there are localised regional signatures among MP industries that may reflect cultural traditions associated with particular populations.

PLOS ONE
The sediment characteristics of this GH suggest a rather warm climate (comparing to the other GHs) leading to a more diverse fauna to appear in the environment. Unfortunately, due to the lack of high-resolution, well-dated paleoclimate records from the Zagros Mountains for the past 45,000 years, correlating our GHs to past regional climate changes is difficult. However, pollen evidence from Lake Zeribar characterised by high percentages of Artemisa and Chenoppodiaceae indicates the dominance of a cool, dry steppe for the period between 40,000 and 14,000 cal. BP. Nevertheless, a slight increase in tree pollen at around 40,000 BP suggests a fairly warmer and wetter (less cold) climate at this time [58], which may correspond to the upper part of GH5.
Sediment compositions of GH4 and GH3 suggesting a prominent cold period characterized by higher snowfall simultaneous with Heinrich Stadial 4 (Figs 3B and 3C and 9). The highly calcareous nature of the sediments might point to an increase in the rate of chemical weathering due to the contact of snow meltwater with the surrounding limestone. One might speculate that the discharge of meltwater runoff rich in calcium carbonate into the rockshelter deposits has caused the highly calcareous sediments here. Moreover, the observed limestone fragments might originate from physical glacial erosion and permafrost degradation. Several lines of evidence such as a deep depression of the snowline to an elevation of about 1800 m in the Zagros Mountains during glacial periods [59], the enhanced rate of calcium carbonate deposition in the nearby Hashilan Wetland during cold climate intervals [60], and extremely sensitive climate of the West Asia to Heinrich cold events [61] provide supporting evidence for our interpretation that the observed highly calcareous sediments in this GH are associated with cold climatic conditions during Heinrich Stadial 4. Pollen evidence from Lake Zeribar reflects an almost treeless environment for the Zagros Mountains between ca. 33,000 and 15,000, which has been attributed to a colder and drier climate [58]. Given a considerable time gap at the transition from GH2 to GH1, it can be postulated that an intervening horizon has most likely been disappeared due to human interferences in the rockshelter environment.
The MP-UP transition in Europe which started between 47 kya BP overlapped with the spread of Homo sapiens and Neanderthals around 40 kya. Recent Palaeogenetic studies show that the gene flow between Neanderthals and Homo sapiens have occurred in roughly between 60-50 probably in southwestern Asia [62,63]. This new information is consistent with the data from Bawa Yawan. The Bawa Yawan Neanderthal suggests an extirpation date of around 45,000-40,000 cal. BP, but the lithic evidence shows extend for a while. If the conditions under which Neanderthals became extinct, locally and overall, are to be determined, it is necessary to develop a more integrated approach to the problem.

Conclusion and implications
The discovery of an in situ Neanderthal tooth at the Bawa Yawan rockshelter in the West-Central Zagros Mountains, its direct association with typical Zagros Mousterian artifacts, and the chronological attempt of the local MP to UP transition, has implications for our understanding of the spatial patterning of Neanderthal extinction and its relationship to expanding Homo sapiens populations and further work need to be done in order to shed more light on this intriguing region. In addition to expose a sequence of three archaeological techno-complexes of MP, UP and Epipalaeolithic formed in four GHs, excavations in Bawa Yawan rockshelter offer evidence that Neanderthals had a demographically dynamic history, shaped by multiple dispersal events and localised population extinctions. Since there is reliable evidence for diverse Homo sapiens in the Middle East during past 45,000, the last Zagros Neanderthals potentially overlapped with groups in space and time.